arXiv:astro-ph/0512357vl 14 Dec 2005 


An investigation of the submillimeter background radiation using 

SCUBA and Spitzer 


S. Dye 1 , S. A. Eales 1 , M. L. N. Ashby 2 , J. -S. Huang 2 , T. M. A. Webb 3 , P. Barmby 2 , S. 
Lilly 4 , M. Brodwin 5 , H. McCracken 6 , E. Egarni 7 , G. G. Fazio 2 


ABSTRACT 

We investigate the redshift dependence of the contribution to the extragalactic far- 
infrared/sub-millimeter background from galaxies detected by the Spitzer Space Telescope at 
8pm and 24pm. Using seven-band optical to mid-infrared photometry, we estimate photometric 
redshifts for the Spitzer sources which appear to be mostly L* galaxies at a median redshift of 
z = 1.0. These sources, extracted from deep 8/im and 24pm mosaics of the CUDSS 14-hour field 
with 5 ct limits of 5.8/iJy and 70pJy respectively, exhibit significant 850pm and 450pm emission 
as observed by SCUBA. At 850pm, after removing > 4 ct sources and those securely identified in 
our previous cross-matching paper, we measure stacked flux at the significance level of 4.4 ct and 
2.9cr from the full 8pm and 24pm galaxy catalogue respectively. At 450pm, flux is detected from 
all 8pm galaxies at the level of 3.5 ct, while there is no significant emission from the 24pm galaxies. 
We find that the 850pm flux is emitted almost exclusively at z > 1.3 from the Spitzer sources 
with 0.44mJy (4.7 ct) per 8pm source and 0.51mJy (2.8 ct) per 24pm source. This corresponds to 
a contribution of (16 ± 3)% toward the 850pm extra-galactic background from the 8pm sources 
and (5.0 ± 1.8)% from the 24pm sources. At 450pm, only the 8pm sources within the redshift 
interval 1 < z < 2 exhibit significant emission with an average flux per source of 3.35mJy (3 .Oct). 
This is a contribution of (37 ± 12)% to the 450pm background. 

Subject headings: infrared: galaxies 


1. Introduction 

Steady advances toward a thorough under¬ 
standing of the population of high redshift sub¬ 
millimeter (submm) sources uncovered by the 
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Sub-millimeter Common User Bolometric Ar¬ 
ray (SCUBA) and the Max-Planck Millimeter 
Bolometer (MAMBO) have been made since their 
detection in the first deep submm survey by Smail, 
Ivison & Blain (1997). The two most controversial 
issues concerning their nature have been, firstly, 
identification of their dust-cloaked energy source 
and secondly, how they relate to local systems. 

The question regarding the energy source has 
now been largely satisfied. The lack of strong X- 
ray emission from these sources (eg. Ivison et al. 
2000b; Fabian et al. 2000; Alexander et al. 2003; 
Almaini et al. 2003; Waskett et al. 2003; Alexan¬ 
der et al. 2005) suggests that their submm flux 
is dominated by re-radiated emission from intense 
star formation, rather than AGN output. Obser¬ 
vations constrain the net contribution of AGN ac¬ 
tivity to the level of ~ 30% (see Chapman et al. 
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2005, and references contained therein). However, 
the question regarding their relationship with the 
local galaxy population remains only partially an¬ 
swered. 

Finding the connection between distant submm 
sources and local galaxies requires an understand¬ 
ing of their evolution. Sources detected in sur¬ 
veys with SCUBA and MAMBO have extremely 
high estimated star formation rates (100 - 1000 
Mq yr _1 ), enough to rapidly form a massive ellip¬ 
tical (Smail, Ivison & Blain 1997; Hughes et al. 
1998; Lilly et al. 1999; Dunne, Eales & Edmunds 

2003) . However, it is not clear from studies of 
the spatial density or clustering of submm sources 
whether a direct link between local ellipticals and 
the submm population can be made quite so con¬ 
fidently (Fox et al. 2002; Scott et al. 2002). 

The main observational limitation in under¬ 
standing the SCUBA population is the relatively 
large beam size. This gives rise to a correspond¬ 
ingly large uncertainty on extracted source posi¬ 
tions and also means that confusion is a concern 
in deeper observations. One of the key ingredients 
needed to help resolve the evolutionary relation¬ 
ship between local and distant submm galaxies, 
determination of submm galaxy redshifts, is there¬ 
fore especially difficult. A successful technique 
is to observe the sources using radio interferom¬ 
etry, as the surface density of radio sources is suf¬ 
ficiently low to make confident associations with 
coincident submm sources. The high astromet¬ 
ric precision in these radio data greatly eases the 
identification of optical counterparts for follow-up 
spectroscopy. Around two thirds of the submm 
sources detected at 850pnr appear to have radio 
counterparts (eg. Ivison et al. 2002; Borys et al. 

2004) , but because the ratio of radio to submm 
flux decreases with increasing redshift (Carilli & 
Yun 1999), radio observations are biased toward 
lower redshift objects (z < 3). This method was 
recently applied by Chapman et al. (2005) using 
Keck spectroscopy to determine the redshifts of 
73 submm sources cross-identified in VLA maps. 
They found a median redshift of 2 ~ 2.2. 

A frequently ignored issue is that the sources 
in the bright SCUBA samples may not be repre¬ 
sentative of the extragalactic far-infrared/submm 
background radiation as a whole. Determining the 
nature of the individual sources that constitute 
the background is of the greatest importance be¬ 


cause approximately 50% of the entire extragalac¬ 
tic background radiation (minus the CMB) is radi¬ 
ated in the submillimeter waveband (Fixsen et al. 
1998). Yet the bright SCUBA sources for which 
there are redshifts may represent only a tiny frac¬ 
tion of this background. The problem is that the 
far-IR/submillimeter background is much brighter 
(~ 30 times greater in terms of ul „) where it peaks 
at ~200pm than it is at 850pm, the wavelength 
of the SCUBA surveys. At 850/im, around 30% 
of the background can be resolved into sources 
brighter than 3mJy (eg. Hughes et al. 1998; 
Eales et al. 2000; Smail et al. 2002; Webb et 
al. 2003b; Chapman et al. 2005). Estimates of 
this fraction extend up to 60% for sources down 
to lrnJy (eg. Smail et al. 2002; Chapman et 
al. 2005). The sources for which Chapman et 
al. (2005) measured redshifts are radio-detected 
sources with 850/im fluxes brighter than ~3 mJy. 
Using estimates of the spectral energy distribu¬ 
tions of individual galaxies Chapman et al. (2005) 
concluded that their sample represents 30% of the 
background at 850pm but only about 2% of the 
background at 200pnr. 

Therefore, despite the success of the work on 
the brighter SCUBA samples, there are many 
unanswered questions. For example, is the redshift 
distribution that Chapman et al. (2005) measure 
representative of the far-IR/submm background as 
a whole? Until the background at 200pm can be 
directly resolved into individual sources, the next 
best alternative is to study the SCUBA sources at 
fainter 850pnr fluxes. Unfortunately, here one hits 
the obstacle of confusion. 

There are two ways round this obstacle. One 
is to exploit gravitational lensing to mitigate the 
problem (eg. Smail et al. 2002). The second 
is the approach we adopt in this paper, using 
the SCUBA images to estimate the submillimeter 
fluxes of objects selected in other wavebands for 
which there are already accurate positions. Pea¬ 
cock et al. (2000) used this reverse procedure by 
stacking 850pm SCUBA detected emission from 
optical sources in the Hubble Deep Field North. 
This study identified significant submm emission 
from galaxies with high UV star formation rates. 
Serjeant et al. (2003) repeated this work replac¬ 
ing the optical data with 15pm ISO sources, but 
measured no emission. However, the much im¬ 
proved sensitivity and resolution of the Spitzer 
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Space Telescope (Werner et al. 2004) recently 
enabled Serjeant et al. (2004) to statistically de¬ 
tect 5.8/tm and 8/tm sources in 850/tm and 450/tm 
SCUBA data. Finally, Knudsen et al. (2005) re¬ 
cently detected significant 850/tm emission from 
distant red galaxies selected by J—K> 2.3, finding 
that these sources are probably strong contribu¬ 
tors to the far-IR/submm background. 

In this paper, we stack 450/tm and 850/tm flux 
observed by SCUBA at the position of 8/tm and 
24/tm Spitzer sources in the Canada-United King¬ 
dom Deep Sub-millimeter Survey (CUDSS) 14- 
hour field. From these data, we estimate the frac¬ 
tional contribution from Spitzer sources to the ex- 
tragalactic background at 850/tm and 450/im. Us¬ 
ing a combination of ground-based optical, near 
infrared and Spitzer 3.6/tm and 4.5/tm observa¬ 
tions, we compute photometric redshifts for the 
Spitzer sources, to investigate the epoch at which 
their attributed submrn flux is emitted. This pa¬ 
per follows on from our previous paper (Ashby et 
al. 2005) where the forward approach of cross- 
identifying Spitzer sources with submrn sources 
was carried out. 

The layout of this paper is as follows: In Section 
2 we describe the data. Section 3 details our pho¬ 
tometric redshift estimation. Stacking is carried 
out in Section 4. We conclude with a summary 
and discussion in Section 5. 

2. Data 

Observations of the CUDSS 14 hour field anal¬ 
ysed in this paper comprise three distinct sets: 1) 
Mid-infrared observations using the Spitzer Space 
Telescope, 2) Submrn data observed with SCUBA, 
3) Ground-based optical and near infrared obser¬ 
vations. 

2.1. Spitzer Space Telescope Data 

The Spitzer observations discussed in this pa¬ 
per were obtained as part of the Guaranteed Time 
Observing program number 8 to image the ex¬ 
tended Groth strip, a 2° x 10' area at a ~ 14 h 19 m , 
S ~ 52°48 / (J2000) with the Infra-red Array Cam¬ 
era (IRAC; Fazio et al. 2004) and Multi-band 
Imaging Photometer for Spitzer (MIPS; Rieke et 
al. 2004). Approximately 90% of the CUDSS 14 
hour held falls robustly inside this area, the re¬ 
maining 10% in the south-east corner having poor 


or no coverage. This corner was cropped in all 
data sets presented here (see Figure 1). 

Carrying on from our identification paper 
(Ashby et al. 2005), we have used the 8/tm IRAC 
and 24/tm MIPS source positions for the stack¬ 
ing. The fraction of sources detected in the higher 
frequency IRAC data that are dusty submrn emit¬ 
ters is expected to be lower than in the 8/tm and 
24/tm data. In addition, the sensitivity of these 
shorter wavelength data is much higher than the 
SCUBA maps so that the stacking tends to sam¬ 
ple noise rather than the detectable population of 
faint sources. Indeed, Serjeant et al. (2004) failed 
to detect significant submrn emission from Spitzer 
sources observed in the IRAC 3.6/rm and 4.5/tm 
channels in SCUBA observations of the Hubble 
Deep Field. The converse is true of the longer 
wavelength MIPS data where the number density 
of sources seen is too low to act as a good probe 
of the submrn emission observed in our SCUBA 
maps. The 5<r point source sensitivity of the 8/tm 
and 24/mr data is 5.8/iJy and 70/tJy respectively. 

Although we made no use of the 3.6/im and 
4.5/tm IRAC detected sources in the stacking di¬ 
rectly, photometry from these channels was com¬ 
bined with the ground-based optical data for com¬ 
putation of photometric redshifts in Section 3. 
The 5 <t point source sensitivity of both the 3.6/tm 
and 4.5/tm data is 0.9/tJy. 

There are 553 8/tm sources that fall within the 
850/tm map area and 511 that fall within the 
450/tm map. Of the 24/tm sources, 141 fall within 
the 850/tm map and 133 within the 450/tm map. 
All of the 24/tm sources are detected at 8/tm. 

2.2. Sub-mm data 

The SCUBA observations amount to 63 hours 
worth of data taken on 20 different nights over 
the period from March 1998 to May 1999 at the 
James Clark Maxwell Telescope (JCMT). The fi¬ 
nal 450/tm and 850/tm maps of the ~ 7' x 6' sur¬ 
vey region were created by combining many ’sur¬ 
vey units’. Each survey unit was observed for ap¬ 
proximately one hour with a 64-point jiggle pat¬ 
tern (ensuring they are fully sampled), nodding 
JCMT’s secondary mirror and chopping by 30" in 
right ascension. 

The data reduction was carried out with the 
SURF package (Jenness 1997) and a noise map 
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created using a Monte Carlo technique to simulate 
the noise properties of each of SCUBA’s bolome¬ 
ters. We refer the reader to Eales et al. (2000) 
for more specific details of the data reduction. 

We work with beam-convolved signal-to-noise 
(S/N) maps in this paper, to maximise source de¬ 
tection. Since the goal is to detect faint SCUBA 
sources buried in the map noise, prior to convolv¬ 
ing with the beam, we subtracted SCUBA sources 
based on their submm S/N and their identifica¬ 
tions in our previous paper (Ashby et al. 2005). 
For the 850/im data, we made three maps with dif¬ 
ferent combinations of SCUBA sources removed to 
assess whether the stacking is influenced by this 
removal process. These three combinations are: 

1) Removal of only the sources listed in Ashby et 
al. (2005) with a secure Spitzer identification. A 
‘secure’ source is defined as having a bright 8/(m 
positionally coincident (r < 10") Spitzer counter¬ 
part with IRAC-MIPS colors indicative of a high- 
redshift dusty source (see Huang et al. 2004). 
Ashby et al. (2005) made 7 secure identifications. 

2) Removal of the secure sources along with the 
remaining 5 sources detected by SCUBA with sig¬ 
nificance > 4tr. 3) Removal of the 12 sources in 
case 2 along with the remaining 6 sources hav¬ 
ing a ‘possible’ Spitzer identification at 8/im. A 
source identification is classified as ‘possible’ when 
selected from multiple co-incident sources. As 
in case 1, the most likely counterpart is chosen 
based on IRAC-MIPS colors. Removing this last 
combination of 18 sources leaves only 2 of the 23 
SCUBA sources in the 850/im data (three fall in 
the cropped south-eastern corner), both with sig¬ 
nificances of 3.0cr (Webb et al. 2003a). 

Sources were removed from the raw SCUBA 
maps by subtracting the beam profile at the posi¬ 
tion of the source, scaled by the source flux. There 
are no directly detectable significant sources in 
the 450/trn data, hence no sources were removed 
from the 450/im map. Figure 1 shows the beam- 
convolved S/N map at 850/mr with the secure + 
> 4cr sources removed and the (unmodified) map 
at 450/xm. The 8/(m and 24/im Spitzer source po¬ 
sitions within the CUDSS 14h area having 2 > 1.3 
are over-plotted on the 850/im map and those with 
1 < z < 2 are over-plotted on the 450/im map, to 
coincide with the findings of Section 4.2. 


2.3. Optical Near Infrared Data 

In addition to the IRAC matched photome¬ 
try, we also matched to optical and near-infrared 
ground-based data. For the optical photometry, 
we extracted U, B, V & I fluxes from the Canada- 
France Deep Fields survey (CFDF) (McCracken 
et al. 2001) imaged with the 4.0m Blanco Tele¬ 
scope and the Canada-France-Hawaii Telescope 
(CFHT). Like the IRAC photometry, fluxes were 
summed in 3" diameter apertures. The 3(7 point 
source sensitivities in AB mags are 27.71, 26.23, 
25.98 and 25.16 for U, B, V & I respectively. 

For the near-infrared photometry, we matched 
to our own SExtracted (Bertin & Arnouts 1996) 
catalogue of the K-bancl image of Webb et al. 
(2003a), observed with the CFHT using the 
CFHTIR camera. The K-band data reach a depth 
of ATab ~ 23. 

3. Photometric Redshifts 

From the matched aperture photometry in U, 
B, V, I & K as well as the IRAC 3.6/im and 4.5/im 
channels, we obtained photometric redshifts for 
the 8/tm and 24/tm sources using the HyperZ red- 
shift code (Bolzonella et al. 2000). We de¬ 
cided not to use photometry from the two longest- 
wavelength IRAC channels, 5.8/im & 8/(m, be¬ 
cause these data are less sensitive than the two 
shortest-wavelength channels, they are more af¬ 
fected by dust emission features (Lu et al. 2003) 
and the spectral templates are uncertain at these 
wavelengths (specifically, the contribution from 
PAHs). We excluded sources with poor quality 
measurements at 3.6/tm and 4.5/zm, although we 
did not exclude sources with poor measurements 
or non-detections in the longer wavelength bands. 

The local galaxy spectral energy distributions 
(SEDs) packaged with HyperZ are those given by 
Coleman et al. (1980), with an extrapolation into 
the infrared using the results of spectral synthesis 
models. The lack of any empirical basis for these 
templates at wavelengths > 1/rni is clearly unsat¬ 
isfactory since our photometry includes both near- 
infrared and mid-infrared measurements. There¬ 
fore, we constructed our own set of templates in 
the following way. Mannucci et al. (2001) list em¬ 
pirical SEDs extended from 0.1/im to 2.4/tm for 
the Hubble types E, SO, Sa, Sb and Sc. We ex¬ 
tended these out to 6/nn using the average SEDs 
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Fig. 1. - The CUDSS 14 hour field showing IRAC 8^m sources (circles) and MIPS 24/mr sources (crosses) 
overlaid on the SCUBA data. Left: The 850/im beam-convolved map with 2 > 1.3 Spitzer sources overlaid 
(lighter shading indicates a higher subrnm flux). In this map, the ‘secure’ and > 4<r SCUBA sources have 
been removed (see text). Right: The 450/im beam-convolved map with 1 < 2 < 2 Spitzer sources overlaid. 
In both maps, the south-eastern (bottom-left) corner is cropped where there is no Spitzer coverage. 


for disk galaxies and elliptical galaxies listed in 
Table 3 of Lu et al. (2003) using the disk-galaxy 
SED for all Hubble types apart from ellipticals. 
The one disadvantage of these templates is the 
lack of a template for an irregular galaxy, and we 
therefore retained the HyperZ template for an Im 
galaxy, extending this into the mid-infrared us¬ 
ing the average disk-galaxy SED from Lu et al. 
(2003). 

We tested the accuracy of our photometric red- 
shifts by estimating redshifts for the 8/im sources 
in the field that have spectroscopically confirmed 
redshifts. There are 61 galaxies from the Canada- 
France Redshift Survey (CFRS) (Lilly et al. 1995) 
that are detected by Spitzer at 8/jm and an addi¬ 
tional four objects that are possible SCUBA de¬ 
tections with spectroscopic redshifts (Chapman et 
al. 2005). Figure 2 shows the photometric red- 
shift estimates plotted against the spectroscopic 
redshifts. The agreement is fairly good with the 
scatter, defined as i/E[Az/(l + z)] 2 /N, equal to 
0.217 for the CFRS sources and 0.161 for the four 
SCUBA galaxies. We note, however, that there 
are very few galaxies with spectroscopic redshifts 


at z > 1. We are therefore less confident about 
the accuracy of estimates beyond this redshift, al¬ 
though the large bin size used when stacking by 
redshift in Section 4.2 will greatly reduce the im¬ 
pact of large errors. 

Figure 3 shows the I-bancl magnitude versus 
photometric redshift for all the objects in our 8/xrn 
and 24/znr catalogue. The curve shows the pre¬ 
diction for a non-evolving L* galaxy, with the K- 
correction calculated using the SED for an Sbc 
galaxy. We have used the value for L* given by 
Blanton et al. (2001). The diagram shows that 
the galaxies detected by Spitzer at 8/xrn are mostly 
L* galaxies with a fairly small dispersion about 
this luminosity. This small dispersion is additional 
confirmation that our photometric-redshift tech¬ 
nique works reasonably well (approximately 50% 
more scatter is seen if the Spitzer photometry is 
omitted in the redshift analysis). 

4. Stacking analysis 

The stacking process determines the sum of the 
subrnm flux detected at the position of Spitzer 
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Fig. 2.— Comparison of photometric redshifts 
(using U, B, V, I, K, 3.6/rrn & 4.5/zm fluxes) with 
spectroscopic redshifts for the 8^im sources in the 
CUDSS 14h field. Open circles show the 61 8/xrn 
sources with redshifts measured by the CFRS. 
Filled circles are the 4 SCUBA sources with spec¬ 
troscopic redshifts measured by Chapman et al. 
(2005). 

sources and its significance. In this work, we cal¬ 
culate the total weighted flux, given by 

f = Npifivr 2 )/^ 2 ), (i) 

where /, is the subnnn flux taken from the beam- 
convolved map at the position of the ith Spitzer 
source, cq is its error and N is the total number 
of sources stacked. 

We measure the significance of the stacking in 
two separate ways: 

- The first uses the Kolmogorov-Smirnov (KS) 
test to quantify the significance of the dis¬ 
crepancy between the distribution of S/N 
drawn from the submm map at the position 
of the Spitzer sources and the S/N distri¬ 
bution of all pixels in the entire map. If the 
former distribution is positively skewed with 
respect to the latter, then this is attributed 
to the detection of faint SCUBA sources. 

- The second is based on the significance of the 
actual total weighted flux measured. The 



z 



z 


Fig. 3.— I-bancl magnitude versus photomet¬ 
ric redshift for all 8/rrn sources (top) and 24/im 
sources (bottom). Solid line shows prediction for a 
non-evolving L* Sbc galaxy. 

formal error on / is defined as 

a f = N/y/Xtar*. ( 2 ) 

One can then define the significance /? = 
f/crf. However, Serjeant et al. (2003) 
showed using Monte Carlo simulations that 
/3 is not normally distributed. We correct for 
this with our own Monte Carlo simulations 
(see below). 

We determined the distribution of f3 and the KS 
significance by carrying out Monte Carlo simula¬ 
tions for each of the four SCUBA datasets consid¬ 
ered in this paper (the three 850/rrn maps and the 
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450pm map). 100,000 realisations of the Spitzer 
source positions were performed per dataset. We 
found that the distribution of KS significance is 
perfectly Gaussian in every case, whereas the dis¬ 
tribution of (3 is skewed positively by an amount 
depending on the submm dataset. 

To account for the non-Gaussianity in /3, we 
used its distribution function calculated from the 
Monte Carlo analysis to make plots showing the 
relationship between (3 and its equivalent Gaussian 
significance, (3 g . All four SCUBA datasets gave 
linear relationships which were fitted as follows: 
850pm secure (3 g = 0.79 (3 — 0.15; 850pm > 4cr + 
secure (3 g = 0.95/3 — 0.09; 850/mr possible + secure 
(3 g = 1.01/3-0.22; 450pm/3 S = 0.86/3 + 0.26. From 
these relationships, it is apparent that the 850pm 
data with only secure sources removed would quite 
heavily over-predict the stacking significance if 
not corrected. The corrections for the remaining 
SCUBA data are fairly modest. 

In this paper, all total weighted flux signifi¬ 
cances quoted are the corrected version (3 g . Abso¬ 
lute errors quoted on the total weighted flux are 
those given by equation (2). 

Despite the fact that the KS significance ad¬ 
heres to Gaussian statistics, we found that it of¬ 
ten gives false peaks in simulations of alignment 
maps described in the following section, whereas 
the corrected weighted flux significance, (3 gi does 
not. For this reason, we use (3 g to assess data 
alignment. However, in Section 4.2 where we dis¬ 
cuss stacking Spitzer sources selected by redshift, 
we calculate both (3 g and the KS significance for 
comparison. 

4.1. Stacking the full 8pm 24pm cata¬ 
logues 

A crucial assumption in the stacking analysis is 
that the astrometric solutions of the Spitzer and 
SCUBA datasets are not offset from each other. 
We therefore tested the alignment of the SCUBA 
and Spitzer data by calculating /3 g for a range of 
offsets in RA and Dec added to the Spitzer source 
positions. The test plots this as a map with the 
x-axis corresponding to the RA offset and the y- 
axis the Dec offset. If the Spitzer data is well 
aligned with the SCUBA maps and if a significant 
detection is made, a peak in the vicinity of the 
origin should be apparent. This peak should be 


singular. If more than one significant peak is seen, 
then the alignment of both data sets cannot be 
established with confidence and the robustness of 
the detection must also be drawn into question. 

Figure 4 shows the maps of (3 g over the 100" 
range of offsets in RA and Dec spanned, for both 
8pm and 24pm Spitzer sources stacked with the 
450pm data and the 850pm data with secure + 
> 4 a sources removed. The contours start at a 
significance of 2er with intervals of 0.5<r. The more 
smoothly varying features of the 850pnr plots rel¬ 
ative to the 450pnr plots reflects the larger beam 
size at 850pnr. 

Stacking the 8pm sources shows a clear detec¬ 
tion at both 850pnr and 450pnr with peak sig¬ 
nificances of 4.4cr and 3.5(7 respectively. The 
850pm peak occurs at an offset of (0", —2") and 
the 450pm peak at (4",0"). The results from 
stacking with the 24pm catalogue are less clear; 
at 850pm, there are two 2.9(7 peaks at offsets of 
(5",-10") and (-8", 12") with a significance of 
2.5cr at (0", 0"). There is no detection of the 24pm 
sources at 450pnr. 

To determine the significance of these offsets, 
we carried out Monte Carlo simulations of the 
stacking at 850pnr and 450pm. For each reali¬ 
sation, we created a noisy synthetic SCUBA map 
and placed sources with fixed S/N by adding the 
SCUBA beam at random positions. After con¬ 
volving the map with the beam, stacking was car¬ 
ried out, offsetting this randomized catalogue by 
varying amounts to produce a f3 g map, like those 
shown in Figure 4. We carried out 200 realisations 
for each of five different source S/N values. 

Figure 5 shows the mean and scatter of the ra¬ 
dial offset of the peak in (3 g from (0", 0") as a 
function of the significance of the peak. We found 
no discernible difference between the 850pnr and 
450pm results hence this plot serves for both wave¬ 
lengths. From the plot, we conclude that the peak 
of 4.4cr seen in the significance map for the 8pm 
sources and 850pnr data, occurs on average in the 
simulations at a radial offset of r = 2.0" ± 1.5". 
The measured offset of (0", —2") is therefore con¬ 
sistent with there actually being no misalignment 
between the 850pm SCUBA map and 8pm data. 
Similarly, at 450pnr, the offset of (4",0) of the 
3.5cr peak is consistent with the expected offset of 
r = 2.6" ± 1.7" for no misalignment. The impli¬ 
cation that the 450pm and 850pnr maps are well- 
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850,ttm SCUBA data, 8/um Spitzer sources 


450/im SCUBA data, 8/tm Spitzer sources 



Fig. 4. — Variation of the total weighted submm flux significance, j3 g , as a function of the offset of the 
SCUBA maps from the Spitzer source catalogue (i.e. add offset to Spitzer position to get SCUBA position). 
Contours start at 2er level and have intervals of 0.5er. Top row shows results for the 8/im Spitzer sources 
and the bottom row shows the 24/im source results. The 850/im map has secure + > 4cr submm sources 
removed; no submm sources are removed from the 450/im map. 


aligned is reassuring; as an additional check, we 
verified that this is the case from astrometry of 
point source calibrators acquired at both wave¬ 
lengths throughout the CUDSS observing. The 
24/im catalogue has the same astrometry as the 
8/im catalogue since both are tied to the Two Mi¬ 
cron All Sky Survey (2MASS; Cutri et al. 2003), 
hence we conclude that all datasets used in the 
stacking are properly aligned. 

To summarise this section, we can make two 
statements. Firstly, there is significant emission at 
850/im and 450/im from all 8/irn Spitzer sources 
on average. From the 24/im sources, the total 
weighted flux has a significance of 2.5<r at 850/im 
but because the significance map is dual-peaked, 


this is a tentative detection. Secondly, the Spitzer 
data and the SCUBA data are well-aligned. 

In the next section, we split the Spitzer cat¬ 
alogues by redshift to investigate the possibil¬ 
ity of this detection being weakened by fore¬ 
ground/background Spitzer sources not truly as¬ 
sociated with the observed submm emission. 

4.2. Stacking by redshift 

Having assigned photometric redshifts to the 
Spitzer catalogues, we can stack sources separated 
into redshift bins. The bin size must be large 
compared to the typical redshift uncertainty. In 
addition, the number of objects per bin must be 









Fig. 5.— Results of Monte Carlo simulations to 
determine the variation of the mean and scatter of 
peak offset from (0",0") in the significance maps 
as a function of the amplitude of the peak mea¬ 
sured. Results plotted here apply to the 8/trn 
sources and either 850/im or 450/trn data. 

equal so that the significance across bins can be 
fairly compared (i.e. to keep Poisson noise con¬ 
stant per bin). Given the accuracy of our pho¬ 
tometric redshifts discussed in Section 3, we split 
the 8/im sources into 6 redshift bins and the 24/im 
sources into 5. This gives an average bin width of 
Az ~ 0.5. Stacking was then carried out for each 
bin at both 850/im and 450/im. 

The total weighted flux from all sources in a 
given redshift bin (expressed in units of MJy/sr) is 
plotted against the median bin redshift in Figure 
6. The 850/im data show a clear detection with 
both 8/im and 24/im Spitzer sources lying at z> 
1.3. With the 450/im data, only the 8/iin Spitzer 
sources produce a notable detection at z ~ 1.5. 

Choosing our most conservative 850/jm dataset, 
that with all secure and possible sources removed, 
and selecting Spitzer sources with z > 1.3, we 
found that the 8/irri sources are detected at the 
level of 4.7a and the 24/im sources at 2.8a. At 
450/jrn, the 8/iin sources selected by 1 < z < 2 are 
detected with a significance of 3.0a. The 24/im 
Spitzer sources selected by 1 < 2 < 2 are not de¬ 
tected. 


These detections were verified with the KS test. 
We applied the same Spitzer source redshift selec¬ 
tions and chose the 850/im data with secure and 
possible sources removed. Figure 7 shows the re¬ 
sults. In these plots, the S/N distribution of the 
full SCUBA map (the ‘map distribution’) is shown 
as an open histogram and the distribution of S/N 
at the redshift selected Spitzer source positions 
(the ‘source distribution’) as a shaded histogram. 
In all plots, bar that for the 450/im + 24/im data, 
excess emission from faint submm sources causes 
the source histogram to be clearly skewed posi¬ 
tively with respect to the map histogram. The KS 
test quantifies the significances of these skews as 
follows: 850/im emission from the 8/im sources at 
4.5a, 850/im emission from the 24/irn sources at 
3.3a and 450/im emission from the 8/im at 2.3a. 
The 24/irn sources are not detected at 450/im. 

Table 1 summarises the results for all versions 
of the submm data. The weighted flux signifi¬ 
cances agree fairly well with the KS significances, 
although on the whole, they are slightly higher. 
The largest discrepancy occurs in the 850/im map 
with only secure sources removed. In this case, the 
> 4a submm sources and/or the submm identifi¬ 
cations defined as possible by Ashby et al. (2005) 
are detected in the stacking. 

5. Summary and Discussion 

We have measured significant 850/im emission 
from Spitzer 8/im and 24/im sources and 450/im 
emission from the 8/trn sources. The 5a point 
source sensitivities of the Spitzer data are 5.8/iJy 
at 8/tm and 70/tJy at 24/im. By computing photo¬ 
metric redshifts for the Spitzer sources using op¬ 
tical, near infra-red and Spitzer 3.6/iin and 4.5/im 
photometry, we have been able to statistically 
identify the epoch from which this submm emis¬ 
sion dominates. We find that the 850/im flux is 
almost exclusively emitted from Spitzer sources at 
redshifts z > 1.3 up to the highest redshift of z ~ 4 
we measure in our sample. Our most conserva¬ 
tive estimate of this flux (i.e. having subtracted 
secure and possible SCUBA sources) is 0.44mJy 
per 8/nn Spitzer source (4.7a significance) and 
0.51mJy per 24/im Spitzer source (2.8a signifi¬ 
cance). At 450/im, the emission appears to peak 
around z ~ 1.5 with an average stacked flux per 
8/im Spitzer source within 1 < 2 < 2 of 3.35mJy 
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850jim SCUBA data, 8pm Spitzer sources 




A _i_i_i_i_i_i_i_i_i_i_ 

0 0.5 1 1.5 2 2.5 


z 



450pm SCUBA data, 24pm Spitzer sources 



Fig. 6.— Variation of total weighted flux of objects binned by redshift (median bin redshift plotted). Top 
left : The 850/un data stacked with the 8/jm sources, where the solid, dashed and dot-dashed lines indicate 
secure, secure + > 4cr and secure + possible SCUBA source subtraction, respectively. Abscissae are plotted 
slightly staggered for clarity. Top right: The 450/uii data stacked with the 8/rnr sources. Bottom left: 850/un 
data stacked with the 24/tm sources. Bottom right: 450/rnr data stacked with the 24/un sources. 


(3. Oct significance). 

We can estimate the contribution our stacked 
flux makes to the extragalactic far-IR/submm 
background using the spectrum of Fixsen et al. 
(1998). According to their measurements, the 
background flux at 850/un is 0.14MJy/sr. The to¬ 
tal weighted flux we measure from 8/un sources at 
z > 1.3 in our 850/un map with secure and possible 
sources removed is 0.022 ± 0.004MJy/sr, a contri¬ 
bution of (16 ± 3)%. The 24/un sources at z > 1.3 
are responsible for approximately one third of this 
with a fractional contribution of (5.0 ± 1.8)%. At 
450/un, Fixsen et al. (1998) measure a back¬ 
ground flux of 0.47MJy/sr. The weighted total 
flux from the 8/un sources with 1 < 2 < 2 in 
our 450/un data is 0.175 ± 0.058MJy/sr, some 


(37 ± 12)% of the background. 

Our estimate of the 850/un contribution to the 
background from 8/un sources is consistent with 
the fraction of (20 ± 8)% determined by Serjeant 
et al. (2004). Within their errors, Serjeant et al. 
(2004) find that the 450/un flux from 8/mr sources 
is sufficient to account for all of the background 
radiation at this wavelength. This is in contrast to 
our findings that indicate a maximum contribution 
of approximately 50%. 

Including the flux of subnnn sources removed 
from the 850/trn data gives a measure of the com¬ 
bined contribution of directly detected and stacked 
sources to the background at this wavelength. The 
secure and possible subnnn sources have a total 
weighted 850/un flux of (62 ± 5)mJy. Adding this 
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SCUBA data 

j KS Significance 

j Subnnn flux (MJy/sr) 

8/tm sources 

24/tm sources 

8/tm sources 

24/tm sources 

850/tm — secure 

4.8a 

2.6a 

0.032±0.004 (172, 5.5a) 

0.011A0.002 (45, 3.9a) 

850/tm — (secure + > 4a) 

4.7a 

3.4a 

0.024±0.004 (172, 5.1a) 

0.007A0.002 (45, 2.7a) 

850/tm — (secure + poss.) 

4.5a 

3.3a 

0.022±0.004 (172, 4.7a) 

0.007A0.002 (45, 2.8a) 

450/tm 

2.3a 

< la 

0.175A0.056 (170, 3.0a) 

0.035A0.027 (41, 1.4a) 


Table 1: Spitzer 8/tm and 24/tm source detections at 850/tm and 450/tm. Redshift selection 2 ; > 1.3 applies 
to the 850/tm data and 1 < 2: < 2 applies to the 450/tm data. Listed are the significances of the KS tests 
and the total weighted submm flux of all objects within the corresponding redshift selection. Quantities in 
brackets list the number of Spitzer galaxies stacked and the significance of the weighted flux, (3 g (note that 
this is not the ratio flux/error due to its non-Gaussian distribution - see text). 


to the stacked 850/tm flux from the 8/tm sources 
gives a total fractional background contribution 
of (29 ± 3)%. Webb et al. (2003c) estimate 
that Lyman break galaxies (LBGs) contribute ~ 
(20 ± 10)% of the 850/tm background. Even mak¬ 
ing the extreme assumption that there is abso¬ 
lutely no overlap between the LBG and Spitzer 
8/tm population (but see Huang et al. 2005, who 
estimate an overlap of ~ 25%) the total fraction of 
the 850/im background that remains unaccounted 
for is at least 40%. 

The average submm fluxes we measure in the 
stacking correspond to fairly regular systems. To 
demonstrate this, we can use the definition of an 
ultra-luminous infrared galaxy (ULIRG) set by 
Clements, Saunders & McMahon (1999). This 
stipulates that a ULIRG must have a luminos¬ 
ity of at least lO n - 4 L 0 measured at 60/trn by the 
infra-red astronomical satellite (IRAS). With an 
assumed IRAS galaxy SED, the expected 850/im 
and 450/tm flux can be calculated for a ULIRG at 
the median redshift of our catalogue, z = 1.0. Fol¬ 
lowing Clements et al. (2004), we take the coolest 
and warmest IRAS galaxy SEDs from the sample 
of Dunne & Eales (2001) to estimate the range 
of submm emission expected. Using the SED of 
NGC958, which is dominated by cold dust at 20K, 
the expected submm flux from an object with a 
60/im IRAS flux of 10 n - 4 LQ at z — 1.0 is 0.77mJy 
at 850/im and 2.06mJy at 450/im (h 10 o = 0.7, 
Q rn = 0.3, A = 0.7). Similarly, the SED of 
IR1525+36 has a much warmer mix of dust with 
temperatures 26K and 57K in the ratio 15:1, and 
predicts a flux of 0.06mJy at 850/im and 0.20mJy 
at 450/im. The average stacked 850/im and 450/im 
fluxes of 0.44mJy and 3.35mJy we measure from 
the 8/tm Spitzer sources indicates that our stack¬ 


ing is sensitive to borderline ULIRGS. Although 
ULIRGS would be considered extreme systems in 
the local universe, at z ~ 1, they are quite the 
norm (eg. Daddi et al. 2005). 

An interesting result is the peak in 450/tm emis¬ 
sion we measure from the 8/tm sources around 
z = 1.5. At this wavelength, the energy density 
of the background is approximately 25% of the 
maximum at ~ 200/tm (compared to only 3% at 
850/tm). Given that the 8/tm sources make up a 
significant fraction of the 450/tm background, it 
would not be unreasonable to expect that the en¬ 
tire far-IR/submm background is dominated by 
systems around this redshift. This is a slightly 
lower redshift than the median redshift z = 2.2 of 
the population of radio detected submm galaxies 
studied by Chapman et al. (2005). However, the 
Chapman et al. sample is substantially brighter 
(by x40) on average than the submm galaxies 
probed by the stacking in our work. This is consis¬ 
tent with the downsizing scenario in which larger 
galaxies generally form earlier than smaller galax¬ 
ies. However, this result could be heavily influ¬ 
enced by difficult to quantify, varying selection ef¬ 
fects between both studies. Clearly, we must wait 
for observations with a high 200/trn sensitivity to 
provide a complete and definitive answer as to the 
nature of the far-IR/submm background. 
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Fig. 7. Distribution of submillimeter S/N at the positions of the Spitzer sources (shaded histogram, 
left, ordinates) compared with the distribution of S/N over the entire SCUBA map (open histogram, right 
ordinates). For the 850/mr data, secure + > 4cr SCUBA sources are removed and z > 1.3 Spitzer sources are 
selected. No SCUBA sources are removed from the 450/xrn data and 1 < z <2 Spitzer sources are selected. 
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